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Lattice effects on photo-excited states in interacting charge frustrated system are examined. Real time dynamics in 
the interacting spinless fermion model on a triangular lattice coupled to lattice vibration are analyzed by applying the 
exact diagonalization method combined with the classical equation of motion. A photo-induced phase transition from 
the horizontal stripe-type charge order (CO) to the 3-fold CO occurs through a characteristic intermediate time domain. 
By analyzing the time evolution in detail, we find that this characteristic dynamics are seen when the electron and 
lattice sectors are not complementary to each other but show cooperative time evolutions. We also find that there are 
threshold values in the optical fluorescence and electron-lattice coupling for emergence of the photo-induced 3-fold 
CO. The dynamics are distinct from those from the vertical stripe-type CO, in which a monotonic CO melting occurs. 
A scenario of the photo-induced CO phase transition with lattice degree of freedom is presented from a view point of 
charge frustration. 


Emergence of inequivalent electronic charge density and its 
long range order in a solid are termed electronic charge order 
(CO). It is observed in a wide class of solids, e.g. transition- 
metal oxides and organic molecular materials. In these sys¬ 
tems, a rich variety of exotic phenomena are found on the 
verge of CO, such as metal-insulator transition, colossal mag¬ 
netoresistance, superconductivity and so on. 31 Regarding 
the real-space CO pattern, lattice geometry plays an essential 
role. In particular, geometrical frustration effects on the inter¬ 
acting charge systems have been long-standing issues. 4,5 * In 
such charge frustrated systems, phase competition and the ex¬ 
istence of a macroscopic number of degenerated charge con¬ 
figurations are expected to affect not only the stabilization of 
CO but also the related phenomena realized in it. 

In the meantime, there have been recent developments in 
artificially controlling the CO states by external field, using 
different experimental tools. Among such external stimuli, ul¬ 
trafast optical pulse has been recognized as a powerful tool 
of the CO control technique. Not only the optically induced 
CO melting, but also emergence of a hidden CO state and 
other exotic optically induced phenomena have been observed 
experimentally. 6 10 * The CO systems are recognized as suit¬ 
able targets in such ultrafast optical studies, since the CO in¬ 
sulators are more susceptible than insulators owing to other 
mechanisms such as the Mott insulators. In the transient op¬ 
tical responses in CO materials, it is widely accepted that lat¬ 
tice degrees of freedom play a key ingredient. First of all, 
the characteristic time scales for the electron and lattice de¬ 
grees of freedom are different, the former being faster than 
the latter in general. Then, the electron-lattice coupling gives 
an energy flow from the electron sector to the lattice sector. 
Furthermore, by separating the electron-electron and electron- 
lattice interactions, having different characteristic time scales, 
transient optical study may identify the dominant stabilization 
mechanism of the CO at the equilibrium. Recently developed 


time-resolved diffraction and spectroscopy techniques greatly 
promote direct observations of the transient lattice dynamics 
in the electron-lattice coupled CO state. In order to study such 
issues, theoretical calculations for the photo-induced CO sys¬ 
tems in which the lattice degrees of freedom are taken into 
account have been performed so far. 6 - 8 - 11 13 * Combination ef¬ 
fects on the electron-lattice interactions and geometrical frus¬ 
tration should be important, since they both underlie many 
CO materials especially molecular compounds. 14 * Neverthe¬ 
less, these effects on the photo-excited transient processes and 
its mechanism are not fully clarified yet. 

In this Letter, we examine lattice effects on the photo-exited 
CO states in a geometrically frustrated lattice. We analyze 
the interacting spinless fermion model on a triangular lattice 
where local lattice vibrations and the Holstein-type electron- 
lattice interaction are taken into account. This model with¬ 
out the lattice degree of freedom show three kinds of COs in 
the ground state, according to the anisotropy of the Coulomb 
interactions, i.e. the horizontal stripe-type (h-stripe), vertical 
stripe-type (v-stripe) and 3-fold CO states [see Fig. 1(b)]. 
15181 Photo-induced electron dynamics in this model were ex¬ 
amined in our previous paper. 191 It was shown that the photo- 
induced phase change occurs from the h-stripe to 3-fold COs, 
which is attributable to the charge frustration effects in the 
photo-excited states. Here, by adding the electron-lattice in¬ 
teraction, we will see in the following that the photo-induced 
transition from the h-stripe to 3-fold COs occurs through an 
intermediate time domain, where the electron and lattice sec¬ 
tors are not complementary to each other, but show coopera¬ 
tive time evolutions. This is in contrast to the photo-excitation 
in the v-stripe CO, as well as that in the h-stripe CO without 
the lattice degree of freedom. Based on the numerical results, 
we present a scenario for the photo-induced phase transitions 
in charge frustrated electron-lattice coupled systems. 
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The Hamiltonian is given by 

= - X + Y V,j,linj 

<ij ) <ij) 

+ ~y Y^ + pl[) - A Y qiHh (1) 
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where c, (cj) is an annihilation (creation) operator for a spin¬ 
less fermion at site i, and n,- = c? c, is a number opera¬ 
tor. We introduce local lattice vibrations where the displace¬ 
ment and its conjugated momentum at site i are written as 

q, and p,, respectively. The first and second terms repre¬ 
sent the fermion hoppings and the inter-site Coulomb inter¬ 
actions, respectively. As shown in Fig. 1(a), anisotropies in 
the fermion hoppings and the Coulomb interactions are in¬ 
troduced as (f, f), and ( V. , V'), respectively. The lattice vibra¬ 
tion with frequency u>i at is represented in the third term, and 
the Holstein-type electron-lattice coupling is given in the last 
term where A is a coupling constant. 

As in our previous work, 19 ' the optical pump pulse is in¬ 
troduced as the Peierls phase into the transfer integral as 
tij —> where A(r) is the vector potential at time 

r, and Rjj is a relative position vector connecting sites i 
and j, which we take a Gaussian form given by A(r) = 
A p e( tJiUtpY 1 exp[-T 2 /(2r 2 )] cos(w p r) with amplitude A p , 
frequency u> p , a damping factor t p , and a unit vector e. 

Steady and transient properties are calculated in finite size 
clusters where the cluster sizes are taken up to A - 4 x 6 = 
24 sites and the periodic boundary condition is imposed. 
The electronic wave functions I'T(t)) are calculated by the 
Lanczos-based exact diagonalization method 2(1211 as in our 
previous work 19 ’ with a time step 6t. As for the lattice part, the 
Newtonian equation given by g, = -coj at qi+oJi at A in/) is solved 
numerically by using the leap-frog method, r/,(r + dr') = 
qi(r ) + 5t’ pi(j + dr'/2), and pAj + 6t' /2) - p;(r - dr'/2) + 
dr 'F At). We introduce a time step dr', and the Hellmann- 
Feynman force F,-(r) = - (T>\ a tqA T ) + A (n/) T where (• • •) T is an 
expectation value calculated by v F(r)). Validity of the present 
treatment for the lattice degree of freedom is checked by com¬ 
paring the results with the ones calculated by the method pro¬ 
posed in Ref. 22. In the numerical calculations, the number of 
the fermions are fixed to be A/2, the parameters are chosen 
as f = f, V + V' - 12f, and t p = 3/f, and the polarization of 
the pump photon is set parallel to the y direction. We choose 
M - 15, dr = 0.01 /f and dr' = dr x ICC 4 which are sufficient 
to obtain the results with high enough accuracy. 23 ' 

Let us start from the physical properties before photo¬ 
excitation. In Fig. 1(c), we plot the charge order parameters 
defined by 0(k) = A~'| 2,■((«,-- \/2))e' kR ‘\ as a function 
of V for u>iat/t = 0.25 and A/t = 0.8. Finite values of O(k) 
are owing to spatial symmetry breaking due to the lattice 
distortions. The result is very similar to the previous results 
without the lattice degree of freedom. In the limiting regions, 
two prototypical stripe COs are seen, i.e. the h-stripe and v- 
stripe COs characterized by the wave vectors, k h = (0, zr), 
and k\/ = ( 2 n/ V3,0), respectively (see Fig. 1 (b)). In between 
the two COs around the frustration points, V/t ~ V'/t ~ 6 , 
the 3-fold CO characterized by k 3 = (2/r/ V3, 2n/3) appears. 
The phase boundary points are almost the same for the cases 
with or without the lattice, while the amplitudes of the di¬ 
agonal and vertical COs are slightly larger here. The lattice 



Fig. 1 . (Color online) (a) Fermion hoppings and Coulomb interactions on a 
triangular lattice, (b) Schematic charge configurations for the h-stripe, 3-fold 
and v-stripe CO structures. Filled and open circles represent charge rich and 
poor sites, respectively, (c) The V-dependence of the charge order parameters 
in the ground state with keeping V+V' = 12 1. Filled and open arrows indicate 
parameter values chosen in Figs. 2(a) and (b), respectively, (d) The optical 
conductivity spectra before photo-excitation in the h-stripe and v-stripe CO 
phases. Bold arrows represent the energies which are chosen as energies of 
the pump photons. Parameter values for the lattice are (Oiat/t = 0-25 and 
Aft = 0.8. 


distortions follow the charge distribution patterns since in the 
steady state they are just proportional to each other. 

The optical spectra before pumping are shown in Fig. 1(d) 
where the h-stripe and v-stripe COs close to the phase 
boundaries are chosen (V/t = 5.4, and 6.5) . The regu¬ 
lar part of the optical conductivity is given by cr aa (a>) = 
-(Aw)-'Im ( V F 0 | j a (oj - r H + Bo + \q)~'f I'To) where |'F 0 ) and 
Eq are the electronic wave function and energy in the ground- 
state, respectively, j a is a current operator with a Cartesian 
coordinate a(= x, y), and q is an infinitesimal constant. We 
adopt the adiabatic approximation in the calculation of the 
optical spectra, where the lattice coordinates are fixed. As the 
pump photon energies in the time evolutions, we chose the 
energies at which the spectra take their maxima shown by ar¬ 
rows in Fig. 1(d). 

From here, we consider the two cases, the h-stripe case (the 
photo-excitation to the h-stripe CO) and the v-stripe case (that 
to the v-stripe CO). The time dependences of the charge or¬ 
der parameters in the h-stripe and v-stripe cases are shown 
in Fig. 2(a) and (b), respectively. In both cases, the initial 
COs are weakened after the photo-irradiation. As shown in 
Fig. 2(a), in the h-stripe case, not only its melting but also 
emergence of the 3-fold CO are observed, where finally the 
dominant CO parameter is interchanged. Characteristically, 
one can see that there are three time domains; the initial CO 
melts monotonically in t < 15/f, amplitudes of all COs are 
small in 15/f < r < 60/f, and finally, 0(ki) develops and 
tends to be saturated in 60/f < r. This characteristic time evo¬ 
lution will be discussed later in more detail. In clear contrast, 
in the v-stripe case, shown in Fig. 2(b), 0(ky) just monotoni¬ 
cally decreases after the photo-irradiation. These transient be¬ 
haviors are not observed in our previous work for the purely 
electronic model, 19 ' and therefore are attributed to the energy 
relaxation to the lattice system. 

Energy flows in the electron and lattice degrees of freedom 
give a hint to understand contrasting behaviors above. Sev¬ 
eral components of energies in the two cases are shown in 
the insets of Figs. 2(a) and (b). We define (7Y C ), and 
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Fig. 2. (Color online) (a) Time dependences of the charge order parameters 
in the h-stripe and (b) v-stripe cases. Real-space charge distributions before 
pumping (rt = -15) and at rt = 120 are also shown in (a). Insets show time 
dependences of several components of the energy. Orange light and purple 
dark shaded areas represent the time intervals when the pump pulses are in¬ 
troduced and the intermediate time domain, respectively. Parameter values 
are chosen to be V/t = 5.4, cjp/t = 8 and A p = 0.8 in (a) and (b), and 
V/t = 6.5, cop/t = 11.2 and A p = 1 in (c) and (d). Other parameter values are 
0)iat/t = 0.25 and A/t = 0.8. 


{'He- 1 ), respectively, as the expectation values of the sum of 
the first and second terms, the third term, and the fourth term 
of the Hamiltonian in Eq. (1). During the time interval when 
the pump pulse is applied, i.e. r < t,i, increases in (7Y e ) and 
{He- 1 ) indicate photo-excitations of the electron sector. Note 
that the conservation of the total energy in t > guaran¬ 
tees good accuracy in the numerical calculations. In the v- 
stripe case, monotonic changes in all components of energy 
after photo-excitations are consistent with the results of the 
charge order parameters shown in Fig. 2(b). On the other 
hand, in the h-stripe case, at the intermediate time domain 
15/f < r < 60/f, {H e - 1 ) shows a plateau, and ( H e ) and {Hi) 
show some basin-like shapes. This plateau behavior in {H e - 1 ) 
implies desynchronization of the electron and lattice sectors. 
Then at r = 60/f, the decrease in {H e - 1 ) is seen, implying a 
reduction of the discrepancy between the electron and lattice 
sectors owing to the compatible lattice distortion associated 
with the emergence of the 3-fold CO. 

Next, we investigate the fluorescence dependences, which 
clearly demonstrate the lattice effects on the photo-induced 
transition from the h-stripe CO to the 3-fold CO. In Figs. 3(a) 
and (c), respectively, we show the averaged order parame¬ 
ters {O dw (k)) in the h- and v-stripe cases as functions of the 
absorbed photon density. We define 0' dw (k) as the order pa¬ 
rameter averaged in the time domain of 120 < rt < 150. 
We measure the absorbed photon density by introducing n p = 
(£ av - Eq)I{Nlo p ), where E dv is the total energy averaged in 
120 < rt < 150, and Eq is the total energy before pump¬ 
ing. In the v-stripe case (Fig. 3(c)), the order parameter de¬ 
creases monotonically with increasing n p as expected. On the 
other hand, in the h-stripe case, a first order transition like be- 
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Fig. 3. (Color online) (a) Fluorescence dependences of the charge order 
parameters after the photo-excitation in the h-stripe case, and (c) those in the 
v-stripe case as functions of the absorbed photon density, (b) A contour map 
of the charge order parameters averaged in 120 < rt < 150 in the h-stripe 
case, and (d) that in the v-stripe case. Pump photon amplitudes are adopted 
in the range of 0 < A p < 2. Other parameter values in (a) and (c) are the same 
as those in Figs. 2(a) and (b), respectively. 


havior occurs from the h-stripe CO to the 3-fold CO around 
n p ~ 0.03, implying the existence of the threshold fluores¬ 
cence. Phase diagrams in the photo-excited states are pre¬ 
sented as contour maps of O dv (k) in the n p -a{= A 2 /{2a>i at t)) 
planes in Figs. 3(b) and (d). The melting of the v-stripe CO is 
suppressed in regions of small n p , whereas it is rather insen¬ 
sitive to the value of a; the melting is mainly of an electronic 
origin. As seen in Fig. 3(b), the 3-fold CO phase appears in 
between the h-stripe CO and a disordered state in large a. 
There is a threshold electron-lattice coupling for the photo- 
induced 3-fold CO. This implies that the electron-lattice cou¬ 
pling stabilizes both the h-stripe and 3-fold COs, although the 
lattice contribution is larger in the 3-fold CO. 

Here we discuss the vibration frequency dependence of the 
time interval where the intermediate time domain appears in 
the photo-excitation from the h-stripe CO to the 3-fold CO. 
In Fig. 4, the time dependences of O(k) for several ai/ a , are 
presented. It is clearly shown that the intermediate time do¬ 
main shrinks with increasing oi/ at . That is, the time interval is 
determined by the vibration frequency, and the lattice vibra¬ 
tion acts as a driving force of the h-stripe CO to the 3-fold 
CO. The decrease in oj/„, hardly affects the amplitude of the 
induced 0(k 3 ), while it enhances the lattice displacement. 

Such pictures are consistent with the time dependences of 
the local electronic charge density, (n,), and the normalized 
lattice displacement, q, = u>i at qi/A. As shown in Figs. 5(a) 
and (b), (n,) and q, in the h-stripe case show almost similar 
time dependence which is featureless in the intermediate time 
domain (We slightly varied the parameters from the results 
above to see the characteristics clearly. Now the time domain 
is 15/f < r < 70/f). On the other hand, in 0 < r < 15/f 
and r > 70/f in the h-stripe case, and in all time regions af¬ 
ter photo-excitation in the v-stripe case (see Figs. 5(c) and 
(d)), q, shows clear oscillation, while oscillatory behavior in 
(«;) is weak. The difference (n,) - q, plotted by green dot¬ 
ted curves in Figs. 5 show small values in the intermediate 
time domain in the h-stripe case, and is distinguishable from 
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Fig. 4. (Color online) Time dependences of the charge order parameters in 
the h-stripe case for several lattice vibration frequencies. Parameter values 
are chosen to be oJiat/t = 0.1, A/t = 0.5, A p = 0.8 in (a), coiat/t = 0.25, 
A/t = 0.8, A p = 0.8 in (b), and GJiat/t = 0.4, A/t = 1, A p = 1.2 in (c). Other 
parameter values are V/t = 5.4 and o) p /t = 8. 


other time domains. Therefore, the electron and lattice sec¬ 
tors do change cooperatively in the intermediate time domain 
observed in the h-stripe case, while, in others, the electronic 
states are almost settled, and the lattice oscillates coherently 
around a stable position, which is almost determined by the 
photo-excited meta-stable electronic states. 

The obtained characteristic time evolutions are interpreted 
qualitatively by the adiabatic potential schemes shown in 
Figs. 5(e) and (f). In the v-stripe case, the photo-excited state 
is a charge melted metallic-like state. In the adiabatic poten¬ 
tial for the photo-excited state, lattices oscillate around the 
stabilization point. On the other hand, in the h-stripe case, 
the Frank-Condon photo-excited state turns into the 3-fold 
CO state. In the translocation processes between the two CO 
states, the electron and lattice degrees of freedom are not com¬ 
plementary to each other, and the electron-lattice interaction 
energy is high. This discrepancy induces a driving force of 
the phase transition, in which the electron and lattice sec¬ 
tors change cooperatively. After the system is settled down 
into the 3-fold CO, the lattice oscillates around the equilib¬ 
rium position in the new adiabatic potential plane. Both in the 
early and late time domains in the transient processes, a pic¬ 
ture for the lattice vibrations on a adiabatic potential is valid 
and the electronic energy scale is much higher than the lat¬ 
tice energy scale. On the other hand, in the intermediate time 
domain, the two energy scales are close with each other and 
the non-adiabatic processes occur. This small electronic en¬ 
ergy scale is attributable to the frustration effect in the photo- 
excited state in the h-stripe C0. 19) 

We discuss a relation of the present results to our previous 
work without the lattice degree of freedom. 19 ’ In Ref. 19, it 
was found that the dominant charge correlation function is in¬ 
terchange from the h-stripe CO type to the 3-fold CO type 
by the photo-irradiation, while the dominant correlation for 
the v-stripe case is only weakened. The dominant difference 
from the previous calculations is the existence of the interme¬ 
diate time domain in the h-stripe case. This can be attributed 
to two facts: 1) the coupling to the lattice degree of free¬ 
dom induces a time lag in the transient excited states, and 2) 
the electron-lattice interaction generates an energy barrier be¬ 
tween the photo-excited Franck-Condon state and the 3-fold 
CO state. 

A number of theoretical calculations for the photo-induced 
electron and lattice dynamics in the CO systems have been 
performed. Not only the photo-induced transitions from one 
CO to another CO but also transitions from CO to metal were 
shown by analyzing the extended Hubbard-type Hamiltoni¬ 



Fig. 5. (Color online) (a) Time dependences of the charge density ((«,)), 
the normalized lattice displacement (qi), and difference between the two 
((«,) - c/i) at a charge rich site, and (b) those at a charge poor site in the 
h-stripe case. Results in the v-stripe case are shown in (c) and (d). Parameter 
values in (a) and (b) are chosen to be V/t = 5.4, = 0.25, A/t = 1. 

iDp/t = 9.6 and A p = 1.2, and those in (c) and (d) are the same as those in 
Fig. 2(b). (e) Schematic adiabatic potentials in the h-stripe case, and (f) those 
in the v-stripe case. 


ans with the lattice degrees of freedom, modeling the actual 
molecular organic salts. 6,8 ' 13,22) In the present study, we adopt 
a more basic model, i.e., the spinless fermion model coupled 
to a simplistic lattice degree of freedom via just the Holstein 
coupling, and focus, in particular, on combination effects of 
the charge frustration and the electron-lattice coupling on the 
photo-induced transient dynamics. One main finding in the 
present study is the existence of the intermediate time do¬ 
main, appearing in the h-stripe case. It is worth noting that the 
electron and lattice sectors show cooperative time evolutions 
in this time domain despite the fact that the original energy 
scales in the two sectors are much different with each other, 
i.e. oji at < tij, Vij. This is attributable to the frustration effect 
in the photo-excited state, which introduces a small energy 
scale in the electron sector. This interpretation is justified by 
the results that the electron and lattice sectors do not evolve 
cooperatively in the v-stripe case where the frustration effect 
is weak. 
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